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The origin of large-scale magnetic field in the universe is one of the greatest mysteries in modern 
cosmology. We present a new mechanism for generation of large-scale magnetic field, based on the 
power-counting renormalizable theory of gravitation recently proposed by Hofava. Contrary to the 
usual case in general relativity, the U{1) gauge symmetry of a Maxwell action in this theory permits 
terms breaking conformal invariance in the ultraviolet. Moreover, for high frequency modes, the 
anisotropic scaling intrinsic to the theory inevitably makes the sound horizon far outside the Hubble 
horizon. Consequently, non-inflationary cosmic expansion in the early universe naturally generates 
super-horizon quantum fluctuations of the magnetic field. Specializing our consideration to the case 
with the dynamical critical exponent z = 3, we show an explicit set of parameters for which (i) 
the amplitude of generated magnetic field is large enough as a seed for the dynamo mechanism; (ii) 
backreaction to the cosmic expansion is small enough; and (iii) the high-energy dispersion relation 
is consistent with the most recent observational limits from MAGIC and FERMI. 
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I. INTRODUCTION 

Hoi'ava recently proposed a class of power-counting 
renormalizable theories of gravity [1]. The power- 
counting (super-)renormalizability stems from the 
Lifshitz-type anisotropic scaling 
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with the dynamical critical exponent z = 3 (or z > 3). 
Because of this scaling, the theory is often called Hofava- 
Lifshitz gravity. Although renormalizability of matter ac- 
tion, e.g. the standard model action, does not require the 
anisotropic scaling, quantum corrections should generate 
terms leading to the anisotropic scaling with a common 
z for all physical degrees of freedom in the ultraviolet 
(UV). For these reasons, in the present paper we shall 
seriously consider the anisotropic scaling with z > 3 for 
matter degrees of freedom, especially for photon. 

Note that the value of z (> 3) in the UV is a part 
of the definition of a theory, provided that the theory is 
renormalizable. Once z in the UV is fixed then terms 
leading to higher z would not be generated by quantum 
corrections. In this paper, for simplicity, we shall restrict 
our consideration to the simplest case where z in the UV 
is 3. However, in principle the mechanism presented in 
this paper works for any values of z in the UV. (See the 
last paragraph of Sec. lVII for a comment on the case with 
general z > 3.) 

Cosmology based on this theory has been investigated 
by many authors, and number of interesting implications 
have been pointed out [1, H, 0, [E H 13 ■ Among them, 
the one particularly relevant to the present paper is that 
the anisotropic scaling of a physical degree of freedom 
leads to a new mechanism for generation of super- horizon 
quantum fluctuations without inflation 

Needless to say, the driving force behind recent en- 
thusiasm for cosmology based on Hofava-Lifshitz grav- 
ity is the fact that this theory is a new candidate for 



quantum gravity. At this moment, it is not yet clear if 
some version of Hoi^ava-Lifshitz gravity makes sense at 
the quantum level and can be applied to the real world. 
Indeed, the version without the proiectability condition is 
already known to be problematic [SI . On the other hand, 
problems pointed out in the literature are absent if the 
projectability condition is maintained and if the detailed 
balance condition is abandoned Therefore, Hofava- 
Lifshitz gravity with the projectability condition with- 
out the detailed balance condition has a potential to be 
theoretically consistent and phenomenologically viable. 
While there still remain many issues to be addressed in 
the future, we may regard Hoi^ava-Lifshitz theory as a 
candidate for the UV completion of general relativity. 
For this reason, it is interesting and important to inves- 
tigate cosmological implications of the theory. 

In this paper we shall focus on the origin of large-scale 
magnetic field in the universe. The observed magnitude 
of the magnetic field at scales of galaxies and clusters is 
about 1/iGauss. At larger scales, on the other hand, there 
is only an upper limit (< 10~^Gauss) from e.g. obser- 
vation of the cosmic microwave background 9]. Galac- 
tic magnetic field can be amplified by the so-called dy- 
namo mechanism [Toj , but a seed magnetic field must be 
provided by other mechanisms in an earlier epoch since 
the dynamo mechanism does not generate magnetic field 
from nothing. As for scales of clusters, efficient ampli- 
fication mechanism from tiny primordial magnetic fields 
to the observed amplitudes has not been established. In 
this sense there is no consensus about how observed mag- 
netic fields at cluster scales could be related to primordial 
ones. However, it has been suggested that, once galactic 
magnetic field is amplified by the dynamo mechanism, 
those amplified magnetic fields can spread over cluster 
scales through galactic outflows (winds/active galactic 
nuclei ejecta) or some plasma instabilities ^12]. For 
these reasons, in the present paper, we shall restrict our 
consideration to the origin of primordial magnetic seed 
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fields responsible for galactic magnetic fields. 

Various generation mechanisms of the seed magnetic 
field have been proposed so far. Among them, scenarios 
based on cosmic inflation are popular (See, how- 

ever, Refs. [13 [Hill, [13, [11 for other scenarios.) This is 
largely because inflation provides a natural framework in 
which physical wavelengths of quantum fluctuations are 
stretched by rapid expansion to scales beyond the Hub- 
ble horizon. However, this (quasi-)standard paradigm 
based on inflation was recently challenged on the basis of 
strong backreaction and strong coupling [19]. If we sim- 
ply demand that the generated magnetic field does not 
backreacts to the dynamics of inflation significantly and 
that the coupling constant in the theory is not extremely 
large then the amplitude of the primordial seeds cannot 
exceed 10~'^'^G in Mpc scales. Since a large coupling con- 
stant makes any perturbative calculations untrustable, 
it is necessary to take the strong backreaction into ac- 
count or/and to analyze the strongly-coupled quantum 
dynamics non-perturbatively. In some model of infla- 
tionary magnetogenesis, the effect of backreaction was 
investigated and it was claimed that the magnetic field 
generated by inflation significantly backreacts to the cos- 
mic expansion in such a way that generation of magnetic 
field is, after all, highly suppressed. At present it is not 
clear whether the strong backreaction or/and the strong 
coupling really spoils other models of inflationary mag- 
netogenesis or not. While it is certainly worthwhile in- 
vestigating this issue in more details, it is also plausible 
to look for alternative mechanisms. 

It is well known that conformal invariance of the stan- 
dard Maxwell action prevents cosmic expansion (includ- 
ing inflation) from acting as a generation mechanism 
of magnetic field [20]: the Maxwell field in the flat 
Friedmann-Robertson- Walker universe does not feel cos- 
mic expansion and behaves as if it were in flat spacetime. 
Therefore, any generation mechanisms need to include, 
one way or another, effects breaking conformal invari- 
ance. 

Interestingly enough, in Hofava-Lifshitz gravity, the 
U{1) gauge symmetry of a Maxwell action permits terms 
breaking conformal invariance. Actually, among them, 
most important in the UV are those associated with the 
anisotropic scaling ((T]). Therefore, breaking of conformal 
invariance is not only possible but also inevitable in the 
UV regime of Hofava-Lifshitz gravity. Moreover, as al- 
ready stated, super-horizon quantum fluctuations can be 
generated without inflation 4]. The essential reason is 
that the sound horizon for high frequency modes is far 
outside the Hubble horizon if a physical degree of freedom 
exhibits the anisotropic scaling. 

The rest of this paper is organized as follows. In Sec.lTTl 
we describe the action for an electromagnetic field in 
Hofava-Lifshitz theory. In Sec. lIIII we describe our mech- 
anism for generation of magnetic fields and present the 
power spectrum of the magnetic field. The results ob- 
tained by qualitative scaling arguments there will be con- 
firmed by explicit calculations in Appendix. In Sec. lIVI we 



shall investigate the backreaction problem raised in Ref. 
^19i] and confirm that the backreaction is small enough 
for a wide range of parameters. Then we shall estimate 
the order of magnitude of the generated magnetic field 
in Sec. |Vl We shall show an explicit set of parameters 
for which (i) the amplitude of generated magnetic field 
is large enough as a seed for the dynamo mechanism; (ii) 
backreaction to the cosmic expansion is small enough; 
and (iii) the high-energy dispersion relation is consistent 
with the most recent observational limits from MAGIC 
[m and FERMI [H. Finally, Sec. EI is devoted to a 
summary of this paper. 



II. ELECTROMAGNETIC FIELD IN 
HORAVA-LIFSHITZ THEORY 

In Hofava-Lifshitz theory, gravity is described by three 
basic quantities: the lapse function N{t), the shift vec- 
tor N'^{t,x) and the three-dimensional spatial metric 
gij{t,x). We can combine these three to form a 4- 
dimensional metric of the ADM form: 

ds^ = -N^dt^ + g^j{dx' + N'dt){dx^ + N^dt). (2) 

The fundamental symmetry of the theory is invariance 
under the foliation-preserving diffeomorphism: 



t-^t'{t), x~^x'{t,x). 



(3) 



This symmetry, combined with the value of the dynam- 
ical critical exponent z (> 3) in the UV, completely de- 
termines the structure of the gravitational action [J, . 

In this paper we investigate the dynamics of the elec- 
tromagnetic field, i.e. a U{1) gauge field, in Hofava- 
Lifshitz theory. The basic quantities are the scalar poten- 
tial Ao(t,x) and the vector potential Ai{t^x). The (gen- 
eralized) Maxwell action must respect the C/(l) gauge 
symmetry as well as the foliation-preserving diffeomor- 
phism invariance. As in the gravity sector, these symme- 
tries, combined with the value of the dynamical critical 
exponent z (> 3) in the UV, completely determine the 
structure of the action. The (generalized) Maxwell action 
is, thus. 



5* = i / Ny^dtd^x 



x(Fo,-7V'Fz,)-G[B,] 



(4) 



where F^i = doAi — diA^, Fij = diAj—djAi, and G[Bi] is 
a function of the magnetic field Bi and its spatial deriva- 
tives. The magnetic field is defined, as usual, by 



D ^ ,Jl „km TP 
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where eijk is the totally anti-symmetric tensor with 
£i23 = \/g- Restricting our consideration to the case 
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where Ai is a free field, G[Bi] in general has the form 



G[B,] = a,B,B' 



(6) 



where oi, 02 and 03 are constants and Vi is the spatial 
covariant derivative compatible with gij. The highest 
derivative term in G'[i?i] is the square of the (z — 1)- 
th derivative of the magnetic field. The (generalized) 
Maxwell action (|4]) is a special case of the vector field 
action considered in Ref. [3|. 

It is easy to see that the scaling dimension of Aq and 
Ai are (z + l)/2 and (3 — z)/2, respectively. This will be 
important for the estimate of the power-spectrmn of the 
magnetic field in the next section. 

One of most important properties of Hofava-Lifshitz 
theory is that in the UV the theory exhibits the 
anisotropic scaling with the dynamical critical exponent 
2 = 3^. As already stated in introduction, this property 
should be shared with matter fields such as the electro- 
magnetic field. Therefore, the function G[Bi] in the UV 
should be dominated by the z = 3 term 

1 
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where M is a mass scale defined by 03 — 1/M*. From 
the stability of the system in the UV, the sign of this 
term is required to be positive. For lower energy scales, 
relevant deformations, i.e. terms with less number of 
spatial derivatives, become important. 

Before closing this section, let us briefly mention ob- 
servational bounds on AI. The electromagnetic field in 
our model has dispersion relation ^ 
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M4 M2 

where kphys is the physical wavenumber and k is defined 
by 0.2 — k/M^. This leads to the energy-dependent pho- 
ton velocity 




O ( %)9) 



^ Theories with the dynamical critical exponent larger than 3 are 
power-counting super-rcnormalizable and, thus, worthwhile con- 
sidering. In this paper, for simplicity, we consider the case with 
2 = 3 only. 

^ Rigorously speaking, each coefficient on the right hand side is 
subject to logarithmic running under renormalization group flow. 
How they actually run has not yet been investigated in full de- 
tails. (See, however, Ref. [2^ for analysis in a simplifled setup.) 
Among them, the running of the coefficient of fc^^j^^ is required 
to be small. In the present paper, we suppose that this is already 
achieved by tuning various coupling constants. As for other coef- 
ficients, argument in the present paper does not depend on their 
precise values and we can treat them as constant unless they 
change by many orders of magnitude. 



For K = 0(1), the leading correction to the photon ve- 
locity comes from the term proportional to k'^hysl-^^ ■ 
this case the MAGIC Collaboration ^ and the Fermi 
GBM/LAT Collaborations [22| give similar lower bound 
on M: 



(10) 



If |k| <C 1 then the leading correction comes from the 
term proportional to k^f^yg/M^ and the lower bound on 
M is weaker. 



III. GENERATION OF SUPER-HORIZON 
SCALE MAGNETIC FIELD WITHOUT 
INFLATION 

Let us consider the electromagnetic field described 
by the (generalized) Maxwell action ([4]) in the flat 
Friedmann-Robertson- Walker (FRW) spacetime. The 
metric is given by 

ds^ = —dt^ + gijdx^dx^ 
= -dt^ + a'^Sijdx^dx^ 
= a^[-dr]'^ + 5^jdx'dx^], (11) 

where a is the scale factor and rj is the conformal time. 
The Latin indices run over spatial index (i = 1, 2, 3). 

Let us first look at the freeze-out condition by as- 
suming the power-law expansion for the background uni- 
verse Ul- As seen easily from the action, the dispersion 
relation for the vector potential will be significantly mod- 
ified as 



(12) 



where uj is the physical frequency and k is the comoving 
wavenumber. In the UV the dynamical critical exponent 
z is 3 but we shall leave it as a free parameter for a while 
until we need to specify it. Hereafter, for simplicity we 
adopt the unit with M — 1. The fluctuation is expected 
to oscillate (freeze out) if a; ^ (a; <C H). Thus if 



(13) 



is satisfied, vector fields first oscillate and then freeze out 
afterwards. For the power-law expansion, a oc t^, the 
condition of Eq. becomes p > 1/z. Throughout this 
paper we consider cases satisfying the above condition 
since in this case super-horizon quantum fluctuations of 
vector fields can be generated without inflation. 

Now we compute the power spectrum of the magnetic 
field in a qualitative way (See Appendix for quantitative 
analysis). It is easy to guess the scale dependence of it 
by using the scaling dimension of the fields. The vector 
potential Ai has the kinetic term of the form 



i J dr,d^x{d^A,)\ 



(14) 
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Here remember that a canonically normalized scalar field 
has the kinetic term of the form i J dr](Pxa^{dri(j))^ 
Thus, the behavior of Ai = Ai/a should be similar to a 
canonically normalized scalar field. For simplicity, we 
consider a power-law expansion 



(15) 



where we see from the definition of the conformal time 
that t ^ a77^+'?, p = q/{l + q) and q = p/{l—p). For this 
background, the ratio of _ff to w becomes 



H Ha"- 



(16) 



When H Lu, all relevant time scales agree. Thus, at 
the sownrf-horizon crossing H ^ u, the power-spectrum 
of the "canonically normalized" field Ai should follow 
from the scaling dimension of Ai as 



H 



Q,-(3-2)/z ^-{q+l){3^z)/z 



(3-2:)/(9z-9-l) 



(17) 



So far, we have been working in the unit with M = 1. 
Noting that Vb has mass dimension four, we can easily 
recover M as 



(24) 



(See Appendix for explicit confirmation of this result.) 

The correlation length of the generated magnetic fields 
is roughly the sound-horizon size. Since the anisotropic 
scaling in the UV regime makes the sound horizon far 
outside the Hubble horizon, the magnetic field on super- 
horizon scales can be generated. Hence the correlation 
length naturally becomes the cosmological scale. 

We find that time evolution of the power spectrum 
of the super-horizon magnetic field is proportional to 
rf'/a^. Here if a is the de Sitter expansion a oc l/ry, 
then we obtain Vb oc . Thus, the generated magnetic 
field rapidly decays. In this sense inflationary universe 
is not good for the generation of the magnetic field in 
our model. Therefore we will not consider inflationary 
phases. 



IV. ABSENCE OF BACKREACTION PROBLEM 



where we used k^ ~ a^^^yy''^^''^^ which holds for H ^ oj. 

From the kinetic term of Eq. p^ . it is easy to see that 
the super-horizon growing mode behaves as 



(18) 



Thus, at the super-horizon scale, the time evolution of 
V^. is given by 



„2 / rrx 2(l-g)/(g2-g-l) 

oc !L oc r;2(i-«) oc (-) 



LO J 



(19) 

where we used Eq. (fTB]) in the last. For the moment, we 
did not take care of the wavenumber dependence. Re- 
covering the wavenumber dependence, then, we obtain 



V^ 



a k -IT, 



2il-q)/{qz-q~l) 



where 



1 



qz-q-1 



qz-q-1 
Thus, we compute VAi\B-^u 

and the power spectrum of the magnetic field is 



Vb 



(20) 



(21) 



(22) 



(23) 



As recently pointed out in Ref. |19j| , we have to check if 
the generated magnetic fields affect the background uni- 
verse. Let us suppose that the z = 3 regime, where mag- 
netic fields at super-horizon scales are generated, begins 
at ?7 = ?7s and ends at 77 = jyout- Here, ?7out is determined 
by -ff (j^out) = M and we shall take the limit ?7s — > — 00 in 
the end of calculation. Then, the total energy density of 
magnetic fields for rj > 7/out is 



(25) 



where fc^(??out) and fc/(??s) stand for the wavenumbers of 
fiuctuations which freeze out at r/out and rjs, respectively. 
Here, for simplicity, we have assumed that k, = 0(1) or 
\n\ <C 1 so that the z = 2 regime is short or absent. 
Then we see that n = 2{5q — 4)/ (2g — 1) > is necessary 
and sufficient for the finiteness of the integral in the limit 
rjs —00. Next, one wonders if electric fields affect the 
background universe. Since Ve VB/k'^, we realize that 
n — 2 > is necessary and sufficient for the finiteness of 
the total energy density of electric fields. 

Together with the freeze-out condition of Eq. (fT5|) , we 
have the constraint for the power of the expansion rate 
as g > 1 (1/2 < p < 1). 

At 77 = 77out , ~ M'^ and the background energy den- 
sity is around H'^M^^ M^M^^, where Mpi ~ lO^^GeV 
is the Planck scale. Then the backreaction from the gen- 
erated magnetic field will be negligible if M <C Mp\. Im- 
portantly, this is compatible with the lower bound (|10p . 

In summary, the backreaction problem does not appear 
in the cases with g > 1 (1/2 < p < 1) and M <^ Mp\. 
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V. THE EVALUATION OF THE GENERATED 
MAGNETIC FIELD AT EQUAL TIME 

In order to estimate the magnitude of generated mag- 
netic fields, we need to specify the FRW background evo- 
lution in the early universe. For simplicity, we assume 
that an oscillating scalar field dominates the evolution 
of the background FRW universe in the early stage and 
then reheats the universe at 77 = T^rh- Thus, we set q = 2 
for 77 < 77rh- We also suppose that T^out < Vrh < ?7cq, 
where 77 = r]cq corresponds to the matter-radiation equal- 
ity. Then the scale factor behaves as 



2 



Vo J \ Vcq J 



(0 



Vl h <'n <Vcq 
Vcq < V 



(26) 



Let r^cross (?7rh < ??cross < f?eq) bc the couformal time at 
which the fluctuation with the wavenumber k re-enters 
the horizon. Then the spectrum at matter-radiation 
equality (77 = 77oq) is given by 



'PsiVeq) = ('^^^ISfi'j PB(77cross), 
V Oeq / 



(27) 



where 



C^Yi^YvBiVr^). (28) 

V 77rh / V across^ 



Setting z = 3 and q — 2 (thus n = 4), Vsiilih) can be 
computed from Eq. ([24)1 as 



(29) 



^rh 



where a := »7cq/('7rh'7o) so that a = arf for rj < 7/111. 

Since 77out < '7rh, the horizon re-entry occurs in the IR 
regime, where the sound horizon and the Hubble horizon 
agree. This implies that 77cross ^ k^^. Thus, after some 
short calculations we obtain 



(30) 



1/2 

where we have used the relations 77oq = acq 770 and 77i.h — 
0-1/2^-1/4 1/2 

^rh «eq Vo ■ 

Now we can evaluate the order of magnitude of the 
magnetic field at 77 = 77oq as 



B{Vcq) ^ \/VBiVcq) 



ViMpc"^ / v: 



M \ 1/3 



■Ipl 



-7/4 



At galactic scales, the primordial amplitude (|3ip can be 
amplified to the observed amplitude of galactic magnetic 



fields by the dynamo mechanism, following the argument 
given in Ref. [2a |. 

On the other hand, if there is no amplification mech- 
anism then the equal-time amplitude (|3ip would corre- 
spond to about lO^^'^Gauss or lower for scales of IMpc 
or longer at present time. This means that intercluster 
magnetic field predicted by our mechanism is too weak 
to be observed directly or indirectly [26]. As for cluster 
magnetic field, as stated in introduction, there is no con- 
sensus about how observed amplitudes could be related 
to primordial ones. It is, however, possible that the pri- 
mordial magnetic fields generated by our mechanism at 
galactic scales could be amplified by the galactic dynamo 
and then spread over cluster scales [ll|, [ll| . 

Finally we comment on the constraint from the Big 
Bang Nucleosyntheis (BBN). Since the abundance of the 
light elements is observed precisely, BBN gives the upper 
limit on the strength of the magnetic fields. The limits 
on the homogeneous magnetic fields on the BBN-horizon 
size 10~''Mpc) are less than 10~^ Gauss in terms of 
today's values [28|. The magnetic field generated in our 
current model is B ~ 10~^^ Gauss on the BBN-horizon 
scale in terms of today's values. Thus we sec that it is 
consistent with the BBN constraint. 



VI. SUMMARY 

We have presented a new mechanism for generation of 
large-scale magnetic field, based on the power-counting 
renormalizable theory of gravitation recently proposed 
by Hofava. Contrary to the usual case in general rela- 
tivity, the U{1) gauge symmetry of a Maxwell action in 
this theory permits terms breaking conformal invariance 
in the ultraviolet. Moreover, for high frequency modes, 
the anisotropic scaling intrinsic to the theory inevitably 
makes the sound horizon far outside the Hubble horizon. 
Consequently, non-infiationary cosmic expansion in the 
early universe naturally generates super-horizon quan- 
tum fiuctuations of the magnetic field. Specializing our 
consideration to the case with the dynamical critical ex- 
ponent z = 3, we have shown an explicit set of parameters 
for which (i) the amplitude of generated magnetic field 
is large enough as a seed for the dynamo mechanism; (ii) 
backreaction to the cosmic expansion is small enough; 
and (iii) the high-energy dispersion relation is consistent 
with the most recent observational limits from MAGIC 
and FERMI. 

As stated in Sec. HI the value of z (> 3) in the UV 
is a part of the definition of a theory, provided that the 
theory is renormalizable. In the present paper, we have 
restricted our consideration to the simplest case where 
z in the UV is 3. However, in principle the mecha- 
nism presented in this paper works for any values of z 
in the UV. For general z, radiation energy density scales 
as p oc a^'^^^^K This means that a radiation dominated 
epoch of the universe has a power-law expansion a oc 
with p = 2/{S + z). This expansion law satisfies the con- 
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dition (fT3)) . or p > 1/z, if z > 3. For this reason, if 
we consider a version of the Hofava-Lifshitz theory with 
z > 3 then magnetic fields can be generated during a 
radiation dominated epoch. Fm'ther investigation of the 
mechanism with general z is certainly worthwhile. 
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APPENDIX A: QUANTITATIVE ANALYSIS OF 
POWER SPECTRUM 



Here we will carefully compute the power spectrum 
of the generated magnetic field. To do so we will take 
a gauge- fixing as usual(For example, see Ref. [13 )■ 
Without loss of generality, we can choose the gauge of 
= 0. Then we may want to take the transverse gauge 
of diAi = 0. For example, we consider the UV action of 
z = 3. In this case one of the field equations is 



(Al) 



where A := S'^^didj. From this we can easily see that 
diAi = holds for all t if we impose diAi = at an initial 
time. In addition, we can see that the gauge condition is 
consistent with the remaining field equations. Thus we 
adopt the gauge of Aq = diAi = hereafter. 

In the UV limit the action with the critical exponent 
z is approximately given by 



Suv I dr] d'^x 



(d^A,)' 



(-1) 



z+l 



{aMf 



:A,A^A, 



(A2) 



From now on we follow the conventional second quanti- 
zation. First we expand the vector perturbation as 



where ei{k,a) (cr = 1,2) is the orthonormal transverse 
polarization vector and the operators and bt satisfy 



the following commutators 

= (27:)' 6„„,S^'Hk-k'), 



br ,bi 

fe',cr' 



hi ,bi 

k,{T k' 



The vacuum is defined by 

5^^|0)=G for fc. 



(A4) 



(A5) 



The mode functions follows the Klein-Gordon normaliza- 
tion as usual: 

= (A6) 

Then the equation for the mode function becomes 

«r + , ^ A^ur = 0, (A7) 

where the prime stands for the derivative with respect to 
the conformal time. Introducing a-s 



(27r)3 

the equation for xj: becomes 



(aM)2 



Xk = 0- 



(A8) 



(A9) 



Here we define that the physical frequency is a; := 
fcV(a^M^-i). 

Hereafter, we assume that the scale factor has the 
power law, a = arj'^^ with q = p/(l — p) > ^ 1) 
which comes from the condition of Eq. (|13p . Then the 
equation for becomes 

X~ + l3kV-^''-'^'Xk = 0^ (AlO) 
where Pk k'^'^ /{aAI)'^^^'^. The solution is given by 

+C2VvH!^^ (-2^.7^,71/(2.)) ^ (All) 

where ly := —l/2{qz — q — 1) and i?^"' is the j/th-order 
Hankel functions. From the normalization and choosing 
the mode function to be the positive-frequency mode in 
Minkowski spacetime at short- wave length limit, Ci and 
C2 are fixed and then 



u{k,T]) 



(27r)3 



(A12) 
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Now we can calculate the power spectrum of Ai which 
is defined as 



{0\A^j:A^j:,\0) ^ {2nf6,/'\k + k')^VA., (A13) 
where 

A^^= f d^'xe'^-^Aix). (Af4) 



Then the power spectrum for Ai will be 

fc3 



'Pa. 



27r2 

fc3 



|(2^)SH' 



Therefore, the power spectrum of the magnetic field is 



Vb ^ -^Va. 



(A16) 



Now we estimate Vai and Vb in 77 — > cx3 when the 
fluctuations freeze out {H 3> oj). The mode function will 



be approximated by 
u{k, 77) 



(27r)3 V 27r 



(A17) 



Then the power spectrum is obtained as 

Va.\h». = ^'^~^r^' (^/r(-^))2 (z?,)". (AI8) 

Finally ■Ps|ff>t^ becomes 



47r3a4 



(aAf)(^-i)/(«"-«-i)fc"\, (A19) 



where n = 5 — z/{qz — q — 1). This result agrees with 
that of Eq. using the scaling law argument in the 

text. 
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